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Abstract 

Emerging applications involving device-to-device communication among wearable electronics re¬ 
quire Gbps throughput, which can be achieved by utilizing millimeter wave (mmWave) frequency bands. 
When many such communicating devices are indoors in close proximity, like in a train car or airplane 
cabin, interference can be a serious impairment. This paper uses stochastic geometry to analyze the 
performance of mmWave networks with a hnite number of interferers in a hnite network region. Prior 
work considered either lower carrier frequencies with different antenna and channel assumptions, or a 
network with an infinite spatial extent. In this paper, human users not only carry potentially interfering 
devices, but also act to block interfering signals. Using a sequence of simplifying assumptions, accurate 
expressions for coverage and rate are developed that capture the effects of key antenna characteristics 
like directivity and gain, and are a function of the finite area and number of users. The assumptions are 
validated through a combination of analysis and simulation. The main conclusions are that mmWave 
frequencies can provide Gbps throughput even with omni-directional transceiver antennas, and larger, 
more directive antenna arrays give better system performance. 


1. Introduction 

Wearable deviees are positioned to beeome part of everyday life, whether be it in the realm of 
healtheare, the workplaee, or infotainment Q, Q. Mobile wearables open up unique ehallenges 
in terms of power eonsumption, heat dissipation, and networking Q. From a wireless eommuni- 
eations perspeetive, wearable eommunieation networks are the next frontier for deviee-to-deviee 
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Network region A 



(a) Example wearable communication net- (b) A finite network located, for instance, in a train car. Small 

work. The user’s smartphone can act as a circles represent wireless devices and large circles represent block¬ 
coordinating hub for the wearable network. ages. Also shown the concept of a blocking cone Ci, and a blocked 

interferer Xj € Ci. 

Fig. 1. Many users with wearable networks like those shown in (a) will be located in close proximity as in (b), creating mutual 
interference. People block some of the interfering signals. 


(D2D) communication Q . Wearable networks conneet different deviees in and around the human 
body ineluding low-rate deviees like pedometers and high-rate deviees like augmented- or mixed- 
reality glasses. With the availability of newer eommereial produets, it seems feasible that many 
people will soon have multiple wearable deviees [[^, as illustrated in Fig. |l(a) Sueh a wearable 
network around an individual may need to operate effeetively in the presenee of interferenee 
from other users’ wearable networks. This is problematie for applieations that require Gbps 
throughput like virtual reality or augmented displays. The urban train ear will be a particularly bad 
environment with a high density of independent wearable networks located in close proximity Q 
as illustrated in Fig. l(b)[ Understanding the interferenee environment is eritical to understanding 
the aehievable rate and quality-of-experience that ean be supported by wearable communication 
networks as well as the feasible density of sueh networks. 

The millimeter wave (mmWave) band eontains a wide range of carrier frequeneies eapable 
of supporting short range high-rate wireless eonnectivity [[8|. The mm Wave band has several 
desirable features which include large bandwidth, eompatibility with direetional transmissions, 
reasonable isolation, and dense deployability. Standards like Wireless HD Q and IEEE 802. Had 


[ 10| have already made mm Wave-based eommereial produets a reality. Wearable networks might 


use these standards or might use device-to-device operating modes proposed for mmWave-based 
next-generation (5G) eellular systems ini, [jig. Short-range mm Wave eommunieation systems 
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usually focus on high-speed wireless eonneetivity to replace eable eonneetions. However, these 
emerging protoeols have yet to prove their effeetiveness in a highly dense interferenee seenario. 
The tool of stoehastie geometry has been extensively used to study interferenee in large wire¬ 


less networks [13|-[16|. Prior work on mmWave-based networks has also used the results from 


stoehastie geometry to analyze eoverage and rate [17|-[19| while modeling the direetionality of 
antennas and the effeet of bloekages. For analytieal traetability, most work assumes an infinite 
number of mobile deviees spread over an infinite area. These assumptions allow the analytieal 
expressions related to the spatial average of the system performanee to be simplified through 
applieation of Campbell’s theorem p0| . Analysis of the outage probability eonditioned on the 
network geometry in ad hoc networks with a finite spatial extent and number of interferers was 


performed in [211, whieh was extended to the analysis of frequeney-hopping networks in [22|. 


The unique ehannel eharaeteristies and antenna features for mm Wave networks, however, 
were not eonsidered in [ |2T] |, [ |22| |. The mmWave ehannel has been studied for the outdoor 


environment |24| and the signifieant effeet of bloekages on signal propagation is well known 
1. In erowded environments sueh as train ears or airline eabins, human bodies are a main 


and signifieant souree of bloekage of mm Wave frequeneies [26|, [27|. This implies that the very 
same users that wear the interfering transmitters aet to bloek interferenee from other wearable 
networks. 

In this paper, we eharaeterize the performance of mmWave wearable communication networks. 
We foeus on networks operating at mmWave earner frequeneies that are eonfined to a limited 
region and eontain a finite number of interferers while not explieitly modeling the impaet of 
reflections within the finite region or at its boundaries. We develop an approaeh for ealeulating 
eoverage and rate in sueh a network. As mmWave systems are likely to use eompaet antenna 
arrays, we assess the impact of antenna parameters, in partieular the beamwidth and antenna gain. 


on the eoverage and speetral efflcieney of the system. Compared with [[T^-[|T9|, [ |2T| , [ |22| , we 
use the same eomputational approaeh as in [ |2T| , p^ , with assumptions on mmWave propagation, 
antennas, and blockage similar to those in [ [T7| , [ fT9| , though we model people - not buildings - as 
bloekages. Compared with our prior work in [[T] where the interferers were assumed to be at fixed 
loeations, this paper also eonsiders interferers that are randomly loeated. We begin by presenting 
an analysis that leads to elosed-form expressions for the eoverage probability eonditioned on 
the loeation of the interferers and bloekages. Then, through a sequenee of assumptions, we find 
the spatially averaged eoverage and rate when the interferenee and bloekages are drawn from a 
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random point process. The assumptions and analysis are verified through a set of simulations, 
whieh involve the repeated random plaeement of the users aeeording to the modeled point 
proeess. 

The organization of the paper is as follows: We introduee the network topology and signal 
model in Seetion |I^ We deseribe the interferenee model and derive expressions for the signal 


to interferenee plus noise ratio (SINK) distribution and rate eoverage probability in Section III 


In Section we provide numerieal results when the users are loeated at fixed loeations. We 
assert the simplifying assumptions for analyzing wearable networks when the users are loeated at 
random loeations in Seetion |V] and, in Seetion |V^ verify through simulations that the assumptions 
have a negligible effeet on the aeeuraey of the analysis. Finally, we eonelude our work and give 


suggestions for future work in Seetion VII 


II. Network Topology and Signal Model 

Consider a finite network region A with a referenee reeeiver and K potentially interfering 
transmitters. The reference transmitter is assumed to be loeated at an arbitrary but fixed distanee 
Ro from the referenee reeeiver at an azimuth angle (f)Q and elevation 'ipQ. The area of the network 
in the horizontal plane is denoted by |^|, so that the interferer density A = iT/|^|. The interfering 
transmitters and their loeations are denoted by Xj, i = 1, 2, ..., K. We assume the referenee 
reeeiver to be loeated at the origin and represent Xi as a eomplex number Xi = where 

Ri = \Xi\ is the distanee between the transmitter and the reeeiver, and = ZX* is the 
azimuth angle to X^ from the referenee reeeiver. For simplieity, we assume that all the interferers 
are on the same horizontal plane that eontains the referenee reeeiver, though our model eould 
also be easily generalized to handle the 3-D loeations of the transmitters. Further, this assumption 
results in the 2-D blockage model that is elaborated next. 

To model human body bloekages, we assoeiate eaeh user’s body with a eirele of diameter IF, 


as illustrated in Fig. 1(b) These eireles as well as the loeation of their eenters are denoted by Bi. 
Like Xi, Bi is represented as a eomplex number so that Bi = \Bi\e^'^^\ where \Bi\ is the distanee 
between the eenter of the human body bloekage and the reeeiver, and Zi?j is the azimuth 
angle to B^ from the referenee reeeiver. In this bloekage model, a transmitter X* is bloeked if the 
direet path from X, to the referenee reeeiver goes through the eirele assoeiated with any Bj or if 
Xi falls within the diameter-IF eirele assoeiated with any bloekage Bj. The user is assoeiated 
with both a transmitter Xj and a bloekage Bi, and it is possible that transmitter Xj is bloeked 
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by its corresponding blockage Bi. This is called self-blocking, a phenomenon that was studied 


in [271, in the context of 5G mmWave cellular system. If there are no blockages in the path 
from Xi to the reference receiver, then we say that the path is line of sight (LOS); otherwise, we 
say that it is non-LOS (NLOS). We associate different channel parameters with LOS and NLOS 


paths, accounting for different path-loss and fading models inspired by measurements [281, [29 j. 
In this paper, we assume that an interferer i is potentially blocked from the reference receiver 
by Bj,j 7 ^ i. Under this assumption that no signal is self-blocked, the following algorithm is 
used to determine which signals are blocked. 

1) Determine C, the set of all transmitters Xi that have no blockages Bj,j i within a 
distance of IU/2; i.e.. 


f W 

£=lx,:jX,-Bf> — Vj^t 

where [X* — Bf is the distance along the horizontal plane between X^ and Bj. 

2) Sort the blockages from closest to most distant, so that |i?i| < |i? 2 | < ••• < \Bk\ 

3) For each iG{l, 2, ..., K], compute the blocking cones (wedge in ID^Bc^ as 

IL \ _ . f W 


( 1 ) 


Ci = { X E A \ \x\ > \Bi\, £Bi — arcsin 


2\Bi 


< /-X < £Bi -f arcsin 


2IR-I 


( 2 ) 


4) For each i E £, determine if Xi is blocked by checking to see if it lies within any blocking 
cone; i.e., if 


(pi ^ U 

then Xi is blocked. 

An illustration of the blocking cone discussed here is shown in Fig. 

While the antenna gain pattern G{(p, ip) is a complicated function of the azimuth angle (p E 
[—TT, vr] and the elevation angle ip E ) f]> to facilitate analysis, we use the three-dimensional 
sectorized antenna model as shown in Fig. We characterize the antenna array pattern with 
four parameters - the half-power beamwidth 9^') in the azimuth, the half-power beamwidth 
9^^^ in the elevation, antenna gain G within the half-power beamwidths (main-lobe) and gain g 
outside it (side-lobe). We use the subscript t to denote an antenna parameter for a transmitter 
and subscript r for the receiver. For example, the main-lobe gain of the transmitter is Gt and that 


'strictly speaking, a blocking cone is an instance of a truncated cone because it does not extend to the origin. 
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TABLE I 

Antenna parameters of a uniform planar square antenna 


Number of antenna elements 

A 

Half-power beamwidth, 9^^^ = 9^^'> 

73 

vW 

Main-lobe gain G 

A 

Side-lobe gain g 




Fig. 2. Sectorized 3D antenna pattern. 


of the receiver is Gf. Similarly, the side-lobe gains are denoted by and g^. Because the system 
model and the analysis presented in this paper are general, substituting the appropriate values 
for the four parameters G, g, 9^^'^ and 9^^'> into the expressions corresponding to the transmitters 
and the reference receiver enables the rapid evaluation of the SINK distribution. To compare 
performance in terms of directivity and gain based on practical antennas, we assume that a 
uniform planar square array (UFA) with half-wavelength antenna element spacing is used at the 
transmitters and the receiver. The number of antenna elements at the transmitter and receiver 
are denoted by At and Ar, respectively. The antenna gain G{(f),'ip) of a UFA is modeled as 
a sectorized pattern as follows. The half-power beamwidths in the azimuth and the elevation 


are inversely proportional to a/A [30|. The main-lobe gain is taken to be A, which is the 


maximum power gain that can be obtained using A-element antenna array. Note that this is 
an approximation, though it is possible to design antennas to give near-flat response within the 
beamwidth with G oc N. The side-lobe gain is then evaluated so that the following antenna 


equation for constant total radiated power is satisfied p0 | 

r 

) 

(7(0,-0) cos('^)d'^d0 = dvr. 



(4) 


By using Q, we ensure the antenna arrays are passive components. The values for the half¬ 
power beamwidths (which are equal in both the azimuth and elevation for UFA), main-lobe and 
side-lobe gains for an A element (i.e. a/A x a/A) UFA are given in Table When the number 
of antenna elements is one, we say that the UFA is omni-directional and, hence, the main-lobe 
and side-lobe gains are unity. This serves as a reference to compare the impact of antenna gain 
and directivity. As in p^ , p^ , we assume that each interferer is transmitting with its main-lobe 
pointed in a random direction. 
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We assume Nakagami fading for the wireless ehannels so that the power gain hi due to fading 
from Xi to the referenee reeeiver is Gamma distributed. We use rrii to denote the Nakagami 
faetor for the link from Xi to the referenee reeeiver, whieh assumes a value of m\_ for LOS 
and rriN for NLOS 0. The path-loss exponent for X^ is denoted as a*, where a* = a\_ if 
Xi is LOS and = ctn if it is NLOS. There are different ways to define the signal-to-noise ratio 
(SNR) in a system with antenna arrays: with and without the antenna gains. We use cr^ to denote 
the noise power divided by the referenee transmitter power as measured at a referenee distanee 
exeluding the antenna gains. While cx^ is inversely proportional the SNR, we intentionally do 
not inelude the antenna gains into its eomputation, so that our results will naturally eapture the 
SNR enhaneement that aeeompanies the use of larger antenna arrays. The transmit power of 
Xi is denoted as Pi. Eaeh interferer transmits with probability p^, whieh is determined by the 
random-aeeess protoeol and user aetivity and is assumed to be the same for all interferers. 

We assume that the referenee eommunieation link is always LOS. The referenee link undergoes 
Nakagami fading with parameter mo = m\_ and has path-loss exponent = a\_. Of eourse, it 
is possible that the referenee user’s body itself will ereate bloekages on the referenee link in a 
wearable network. When this oeeurs, it ean be handled in our model by setting mg = and 
cto = q;i\|. Capturing self-bloekage of the referenee link in a more refined model and ineorporating 
the results into the analysis is an interesting topie for future work. 

It is important to note that the boundaries of the finite area are assumed to be impenetrable, 
so there is no leakage of external interferenee into the finite area. Further, refleetions due to the 
boundary and objeets within the network are not explieitly ineorporated in the model. They are 
aeeounted for only in a eoarse way in the different LOS and NLOS model parameters, whieh 
ideally would be determined based on ray traeing or measurement results. The assumption of 
omitting refleetions holds true in many seenarios where the boundaries of the finite area are 
made of poorly refleeting materials sueh as eonerete or brieks. 


III. Interference Model 


Conditioned on the network (meaning the loeations of the transmitters and bloekages), we 
ean find the eomplementary eumulative distribution funetion (CCDF) of the SINR (also ealled 
SINR eoverage probability [ 17|) by adapting the analysis in pT| , [22|. The analysis that follows 
in Seetion |III-A is very general sinee it ean admit the individual interferers to have separate 
and independent values for the ehannel parameters - a* and m*, and does not require the FOS 







channel to have values (ol, mC) and the NLOS ehannel to have values (on, mM). The assumption 
of fixing the ehannel parameters of the LOS and NLOS interferers yields traetable analytieal 
expression for spatially averaged SINK eoverage probability in Seetion |V| 

We define a diserete random variable for i = {1,iT} that represents the relative power 
radiated by Xi in the direetion of the referenee reeeiver. With probability (1 —pt), Xi does not 
transmit at all, and henee /* = 0. Otherwise, the relative power will depend on whether or not the 
random orientation of X/s antenna is sueh that the referenee reeeiver is within the main-lobe. 
We assume a uniform orientation of X/s antenna, so that the azimuth angle (j) is uniform in 
[0,27r) and the elevation angle -0 has a probability density funetion (pdf) ^ cos(V’) in f]- 
The pdfs ean be derived by noting that the surfaee area element of a unit sphere is cos (i/j) d0d^, 
a funetion of the elevation angle ^|J. Thus the probability that the referenee reeeiver is within the 


0(3) 


interferer’s main-lobe is -j— sin ( -^ ) = pm- It follows that 


gCs) 


h = 


0 with probability (1 — pt) 
Gt with probability ptPM 


(5) 


Pt with probability pt (1 — Pm) 

Note that a similar approaeh was used in for modeling adjaeent-ehannel interferenee in 
frequeney hopping: when the interferer transmitted, one of two power eompensations was applied 
depending on whether the interferer hopped into the same or an adjaeent ehannel. In the wearable 
network eontext, we ean justify randomizing the orientation angles of the interferers beeause: 

(1) the user itself may be randomly moving the orientation of its deviees while using them, and 

(2) the user may have a wearable network with several deviees with different orientations and 
random aetivity, though we assume the medium aeeess protoeol (MAC) of a user’s wearable 
network allows only one of her deviees to transmit at a time. It may be noted here that this 
kind of wearable network is still under development, so the exaet MAC protoeol has not yet 
been deeided. We make a reasonable assumption that the network of a given user is eoordinated 
sueh that only one deviee transmits at a time, while the deviees of different users are not so 
eoordinated and ean therefore eollide. 

Now, let us define the normalized power gain from Xj to be 


n, = 


Pi G~°‘i 


if < 0. - 00 < ^ 




-OLi 

'i 


otherwise 


( 6 ) 
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where a* = on if is NLOS and a* = ai if Xj is LOS. This is the worst-ease situation when 

(e) (e) 

iV’ol < If iV'ol > ^7 we have Oj = V i whieh is a simpler trivial ease. For the 

rest of the paper, we assume the non-trivial worst ease and all the analysis presented hereafter 
extends easily for the trivial ease. The SINK is 


7 = 


Gthoflo 


K 


(7) 


cr^ + ^ lihifli 


2=1 


where Oq = GrRQ^° is the normalized power gain from the referenee transmitter, as we assume 
the referenee transmitter is always within the main beam of the referenee reeeiver. The effeet 
of misalignment of beam in the referenee link was eonsidered at lower frequeneies (e.g. UHF) 
in prior work [ [3T| , [32|. In the wearable eommunieation network eontext, however, sinee the 
distanee of the referenee link is short relative to the beamwidth of the antenna, pointing errors 
will not seriously degrade performanee. (For instanee, with our seetorized antenna model, the 
beam eould be off by half the beamwidth without ehanging performanee.) 


A. Coverage Probability 

Denoting Cl = [flo, ■■■, ^k], the eoverage probability Pc(/9, for a given is defined as the 
CCDF of the SINK evaluated at a threshold 13 and is given by 


P,{f3,n)=F[y>(3\n]. 


( 8 ) 


Substituting (|7]) into @ and rearranging leads to 


K 


Pc{f3,n) 


P 


S > cr^ 




Cl 


(9) 


i=l 

where S = (3~^GthQCo, and Yj = Jj/iiDj. Conditioned on Cl, let /Y(y) denote the joint pdf of 
(Yi,..., Y/f) and /s(s) denote the pdf of S. Then, Q ean be written as 


PcW,^) 




fs{s)ds 


/Y(y)dy. 


( 10 ) 


Defining (3o = Pmo/GtClo and assuming that mo is a positive integer, the random variable S is 
gamma distributed with pdf given by 


/s(s) 


(mo-l)! 


s > 0. 


( 11 ) 
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Using ([TT]), the inner integral in (fTO]) is 


/ mQ—l , „ 2\e ( 1 ^ 

/sWd. = 1+1^,. 

^=0 ■ V j=l 


( 12 ) 


<^2+Ef i/i 

Substituting ( [T^ into ( fT^ leads to 

^=0 


Fc(Af^) = e- 


i\ 


y--ye fY{y)(iy- 


(13) 


Using the binomial theorem followed by multinomial expansion, 


K 




y^ = 


i=l 


t=0 


K 


E UE 


Vi = 


i=l 


E 

t=o 


A t\ 


t J a 


2t 


E n 

St \i=l 


K U 

yl_ 

tA 


(14) 


where the last summation is over the set St eontaining all length-ii" non-negative integer se- 
quenees {U,..., tx} that sum to t. This ean be pre-eomputed and saved as a matrix as explained 
in [ |2T| . Substituting ( [T?] ) into ( fT3| ) gives 




PcW,n) = e 




£=0 


i\ 


t=0 


A t\ 


t J a 


2t 


E 

-St ■ 


JJ ^g-/3oi/tjY(y)dy. (15) 


2 = 1 


Given f2, the {Yj}^i are independent. So, f^iy) may be written as Hili/ y*(?/*)> where 


__ V TUi rrii — l 


r(mi) 


Pm- 


HVj 


+ (1 — Pm) 


_rnjJH 

e 3t^i 


u{yi) + (1 -Pt)S{yi), (16) 


Gr* ' - ■ Pt" 

5{yi) is the Dirae delta funetion, and u{yi) is the unit step funetion. From the independenee of 
the {Yj}, ( [15] ) may be written as 


PM(l) = 


£=0 


t=0 


K 


t J a 


2t 


I n ) ’ 


St \i=l 


where 


fco ii 

gtX^t) = / ^e-^°y^fyMdyt. 

Jn 


(17) 


(18) 


To evaluate (18), we use the faet that 


6'=r(fc) 


is a probability density funetion (of a gamma- 


distributed random variable Z) with parameters k,b > 0, so that 


^oo 


b>^T{k) 


dz = 1. 


(19) 


Aeeordingly, ( [T8| ) simplifies to 

u A r(mj -f tj) 
tAT{mi) 


Guigti) =pt[ — 


\pMQtXGt) + (1 - Pm) QtXyt)] + (1 -Pt)(5[A- (20) 


















In ([20|), 6[ti] is the function defined as 
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S[ti\ = 


and Qt, (x) = 1 + 



l^oxVti 

rui 


-(rrii+ti) 


( 21 ) 


( 22 ) 


The assumption of an integer value for mo is key to the derivation of the exact expression for 
the SINK coverage probability in ([T7|). When mo is not an integer, such an exact evaluation is 


not possible to the best of our knowledge. Only an upper-bound using the results from [33| can 
be obtained for a general real-valued mo- 


B. Ergodic Spectral Efficiency 

When the SINK is (3, the spectral efficiency in bits per channel use is 

d = log2(l + /?). (23) 

The CCDF of the spectral efficiency is found by defining the equivalent events 

{7 > /3 I fl} {log2(l + y) >T] \ n} . (24) 

__ > 

■v* 

{7>2'J-1 I n} 

The event on the left corresponds to the coverage probability Pc{(3, f2), while the event on the 
right corresponds to the CCDF of the spectral efficiency, fi), also called the rate coverage 
probability for a given fl. Since equivalent, the two events have the same probability, and it 
follows that 




(25) 


See also Lemma 5 of [ [T7| . 

Using the fact that, for a non-negative X, E[X] = — F{x))dx (see (5-33) in [34|), the 

ergodic spectral efficiency conditioned on f2 can be found from 


EM= /rp.(2--i.n)d„= (26) 

where the last step uses the change of variables (3 = 2^ — 1 dp = j^^^d/?/(l -|- (3). 

In practice, there is a maximum and minimum rate, and hence, a maximum and minimum 
SINK thresholds /3max and /3min, respectively. This maximum may be imposed by the modulation 
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TABLE II 

Antenna Parameters 


Number of antenna elements 

1 

4 

16 

Half-power beam width in the elevation and azimuth (in degrees) 

360 

49.6 

24.8 

Main-lobe gain (in dB) 

0 

6 

12 

Side-lobe gain (in dB) 

0 

-0.8839 

-1.1092 


TABLE III 

Parameters used to obtain numerical results for fixed geometry 


Parameter 

Value 

Description 

Ro 

0.3 m 

Reference link length 

00, 00 

0° 

Antenna main-lobe orientation of the reference receiver 

nil 

4 

Nakagami parameter for LOS link 

rriN 

2 

Nakagami parameter for NLOS link 

a\_ 

2 

Path-loss exponent for LOS link 

On 

4 

Path-loss exponent for NLOS link 

W 

0.3 m 

Width of the human-body blockages 

A 

-20 dB 

Noise power normalized by reference transmitter power 

K 

36 

Number of potential interferers 


order of the eonstellation used and distortion limits in the RF front end while minimum due to 
the reeeiver sensitivity. In this ease, the limits of the integral are /3min and /3max, and 

Pc (/S, n) 


<• 13 ^ 


E[r]] = 


log(2)(l + /3) 


d/3. 


(27) 


The quiekest way to eompute p7| ) is to simply eompute Pc (/3,17) for a finely spaeed /3 and 
then use the trapezoidal rule to numerieally solve the integral. 


IV. Numerical Results for fixed geometry 

In this seetion, we provide numerieal results for eoverage probability and ergodie speetral 
effieieney. The users are loeated at fixed loeations, but to enable a eomparison against random 


topologies (see Seetion VI), their plaeement is eonfined to an annulus A having inner radius r\„ 
and outer radius Tout- Conditioned on the fixed loeations of the interferers and the bloekages, 
the exaet expression for the SINR eoverage probability ean be deriverd using ( [T7] ). We assume 
there are K interfering transmitters, negleet self-bloeking, and assume that the bloekage and 
transmitter assoeiated with eaeh user are eo-loeated; i.e., Bi = Xj for eaeh i. It is assumed that 
the Pi are all the same; i.e., all transmitters transmit at the same power. 
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The values of the antenna half-power beamwidths, main-lobe and side-lobe gains are summa¬ 
rized in Table Note that it is possible to get desired side-lobe isolation by earefully designing 


the array response via windowing similar to filter design This would also add complexity 
to the array design and configuration. Since power and heating issues are critical for wearable 
devices, it is yet to be determined if such techniques would indeed be considered in future 
gadgets. To quantify the effect of antenna directivity and as an example, we chose a uniform 
planar array described in Table The network and signal parameters used to obtain the results in 
this section are summarized in Table |II^ The Nakagami parameters and the path-loss exponents 
assumed are the ones used in (n). For simplicity and to ensure the interferers are uniformly 
spread out in the network region, we let the user locations to be on a n x n square lattice 


restricted to the annulus A. The network region under this assumption is shown in Fig. 3(a) 


where the 7x7 grid locations and the user locations in the network region A are shown. We 
let the minimum distance of two nodes in the grid be 2Ro. For example, when the lattice points 


are separated by 0.6 m (Rq = 0.3 m as in Table and n = 1, we get K = 36 with = 0.3 
m and rout = 2.1 m, which corresponds to an interferer density A = 2.25 (passengers/m^), a 
typical density scenariothat approximates the peak-hour passenger load in urban train cars Q. 

Fig. 1^ shows users placed according to Fig. |3(a)| along with the blocking cones (Fig. |3(b)| ) 
assuming that each user is associated with a blockage of width W = 0.3 m. The blocked users 
are indicated by filled circles. We next provide numerical results for this fixed geometry. The 
dependence of coverage probability on the transmission probability pt of the interferers for a 
fixed transmitter and receiver antenna array configuration is shown in Fig. for the case when 
the transmitters and the receiver use omni-directional antenna. It is seen that, as expected, a 
higher value of pt leads to lower coverage probability for a given SINK threshold. We observe 
similar results for other antenna configurations as well. 

The CCDF of spectral efficiency for different antenna configurations is shown in Fig. for 
a given random-access probability. Here we let pt = 1. Clearly, using more antennas at the 
transmitters and the receiver results in significant improvement in the rate. This is because 
larger antenna arrays provide more directed transmission and reception, thus improving the 
SINK due to the increased antenna gains of the reference link as well as the reduced beamwidth 
of the interfering receivers, which reduces the likelihood that the reference receiver falls within 
a randomly oriented reciever’s main-lobe. The ergodic spectral efficiency for various antenna 
configurations when pt = 1 is summarized in Table |IV| It can be seen that having larger is 
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(a) The locations of the users in a uniform grid of size 7x7 restricted to an annulus. The twelve 
users located outside the circle are deleted from the network. 



(b) The blocking cones associated with the blockages and the blocked users (filled circles). 

Fig. 3. The fixed geometry considered in Section m and the blocking cones associated with the users. The reference receiver 
and the projection of the transmitter onto A are shown in blue and green, respectively. 
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Fig. 4. SINR coverage probability when the users are placed in the fixed positions indicated in Fig. |3(a)| for different transmission 
probabilities pt A't = A^r = 1- Larger pt results in smaller coverage probability. 


TABLE IV 

Ergodic spectral efficiency for various antenna configurations 



I 

4 

16 

I 

0.1762 

0.8710 

1.5481 

4 

1.0880 

2.3282 

3.2820 

16 

2.6734 

4.2190 

5.2850 


more advantageous than having a larger for the fixed geometry eonsidered in this seetion. 
This ean also be seen from Fig. whieh is a eontour plot of the ergodie speetral effieieney 
as a funetion of iVt > 4 and > A for a random-aeeess probability pt = 0.5 and found by 
interpolating the eomputed values of the ergodie speetral effieieney for all integer values of 
and Nr from 4 to 16. We attribute the asymmetrie behavior with respeet to and in Fig. to 
the faet that the interferers have smaller probability of pointing their main-lobes to the referenee 
reeeiver when is large. 

V. Spatial Averaging for Random Geometries 

Infinite-sized networks are usually analyzed by assuming the interferers are drawn from a point 
proeess, then determining the eoverage and rate of a typieal user by averaging over the network 
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Fig. 5. CCDF of spectral efficiency for different antenna configurations when pt = 1 and the fixed network geometry in Fig. 
|3(a)[ Spectral efficiency is improved significantly with more antennas. 



Fig. 6. Contour plot of ergodic spectral efficiency as a function of Nt and Nr when pt — 0.5 for the fixed interferer placement 
according to Fig. |3(a)| 


geometry. Intuitively, this ean be thought of as the performanee seen by a user that wanders 
throughout the network, and thus sees many different network topologies. In this seetion, we 
outline the approaeh we take to analyze the interferenee in a finite sized mm Wave-based deviee- 
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to-device network when the users are loeated at random loeations. While in reality, users are 
generally spaeed far enough apart that their bodies don’t overlap, for mathematieal traetability 
we assume that the users are independently plaeed (whieh inelude eases with overlaps). 

The spatially averaged CCDF of the SINK ean be derived by taking the expeetation of Pc(/^, 
(( fTT] ) in Seetion |I^ with respeet to Q, 


Pc(/5)= P[7>/9]= En[Pc(Af^)]. 


(28) 


This can be performed in two ways: (1) through the use of simulation, and (2) analytically. 
In the former method, we can determine the coverage and rate as follows. We randomly place 
the K potentially interfering users drawn from a binomial point process (BPP) and compute the 
corresponding coverage and rate for each network realization. This is repeated a large number of 
times to obtain the spatial average. While in the limit of an infinite number of trials, this approach 
provides the exact spatially averaged performance, the downside is that it is computationally 


expensive to obtain. The second method is similar to the approach in [^, which however only 
considered Rayleigh fading for the links. The spatially averaged outage probability is found in 


this approach by unconditioning the results we obtained in Section III that were conditioned on 
the location of the interferers and the blockages. Using this approach, we develop closed-form 
expression for the spatially averaged CCDF of the SINK, which is then validated against the 
results obtained via simulation. 


A. Assumptions 

Taking the expectation is complicated by a number of factors that arise primarily due to 
the coupling of interferers and blockages, since each user is not only a potential source of 
interference, but is also a potential source of blockage. To make the problem more tractable, 
we adopt a sequence of assumptions, with each assumption building upon the previous one. 
Simulation results show the validity of the assumptions. 

Assumption 1: The locations of the blockages and interferers are related hy an orbital 
model. Even if a user Bi (which also denotes blockage) is in a fixed location, its transmitter Xi 
could be randomly positioned around it. Hence, we specify the location of in the 2-D plane 
relative to Bi by placing it randomly on a circle of radius d > W/2 and center Bt. Self-blocking 
is now inherent in the model, i.e., if X* is behind Bi, then it is blocked. We refer to this model 
as the orbital model. This is illustrated in Fig. 
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'X,: 



radius-d circle 

Fig. 7. An illustration of the orbital model with the locations of blockage Bi and the interferer Xi 


Assumption 2: The locations of the blockages and interferers are drawn from independent 
point processes. Though this assumption is not by itself that useful, it is a stepping stone towards 
a traetable analysis. With this assumption we can still obtain the corresponding coverage and 
rate using the aforementioned simulation, only now the simulation can lay out K interferers and 
K blockers independently. 

Assumption 3: The blockage states of the interferers are independent. This assumes that 
there is no correlation in the blockage process, even though in reality a transmitter that is close 
to a blocked transmitter is likely to also be blocked. With this assumption we first determine 
the blockage probability Ph{r), which gives the probability of blockage arising from other users 
as a function of distance r to the reference receiver. We can determine Pb{r) either empirically 
(through running simulations of the blockage process) or by using results from random shape 


theory [25|. Then, having established Ph{r), we run a new simulation whereby we first place 
the interferers according to a BPP, then we independently mark each interferer as being blocked 
with probability ph{r). Note that now, we need not place the blockages in the simulation. 

Assumption 4: All interferers beyond some distance Rb are NLOS and those closer than 
Rb are LOS. Here, we replace the irregular and random LOS boundary with an equivalent 
ball. The value of Rb can be found by matching the first moments (Criterion 1 of p^), or 
alternatively, it could be found by matching the average rate. 


B. Analysis of Blocking Probability 

Lemma 1: When the network region A is an annulus with inner radius r\„ and outer radius 
Tout and blockages have diameter W, the probability that an interferer at distance r from the 
reference receiver is blocked by any of the K blockages that are independently and uniformly 
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distributed over A is 

Pb(r) = 

where 


1 - 1 - 




1^1 


1 I 1 rW—^i+u 

^ ^ |yl| 


if Tin < r < Tout - ^ 


K 


(29) 


if r-out - ^ <r < Tout 


p = 


V = 



lyV 2 

— I + Tin arcsm 


W 

2rin 


(30) 




arcsm 




rW 


Tout arccos 


2rr 


out 


r, / ^^OUtN Tout + ?" + W^/2 

- 2y s(s - r)(s - y )(s - —); s =-^-. 


(31) 


Proof 1: See Appendix |A| 

The distance-dependent blockage probability Pb{r) is shown for an annulus of inner radius r\n = 1 
and outer radius rout = 7 in Fig. along with that obtained via simulation. The simulation results 
closely match the analytical approximations derived using Assumptions 1-3. Also, as expected, 
the probability that an interferer is blocked increases as the distance r between the interferer and 
the reference receiver is increased. The dependence of ph{r) on K, the number of interferers 
is shown when IF = 1 in Fig. |8(a) and the dependence on the width of the blockages W is 
shown when iT = 36 in Fig. |8(b) It is seen that with larger values for K and W the blockage 
probability is more for any given distance r from the reference receiver. 

To evaluate Rq under Assumption 4, we compute the mean number of interferers that are not 
blocked as 


K 


^T’out 


p = 27r-^ / {1 — p\^{r))rdr. 


1^1 


(32) 




Equating this average number of non-blocked interferers to the number of interferers in an 
equivalent LOS ball of radius Rq leads to the expression 

/•Tout 


Rb = 


(1 - Ph{r)) rdr + r\„‘ 




0.5 


(33) 


Now Pb{r) is approximated by p\^{r), which is a step function with a step up at distance Rq, 

0 if rin < r < i?B 
1 if < ^ < Tout 


Ph{r) = 


(34) 
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n 

a. 



(a) Pb{r) vs r for different values of K when W = 1. 



r 


(b) Pb{r) vs r for different values of W when K = 36. 


Fig. 8. The distance-dependent blockage probability as K and W are varied. Here we assume = 1 and Tout = 7. 
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C. Analysis of Coverage Probability 

Under assumptions 1 - 4, we ean derive the probability distribution of f^fw) of U* whieh 
are now independent random variables that depend on the loeation of Xi for {1, 2,..., K} 
follows 

(35) 




where 


and 


01 (0*) = 


Oi if|0i-0o|<V 


otherwise 


(36) 


Q2{Ri) = 


a\_ if Ri < 

tt|\i if Ri > Rq 


Oi = ;^G'r and bi = Additionally, we have 


(37) 


(38) 


m-L if Ri < Rb 

rrii = j 

rrtN if Ri > Rb 

For Xi drawn from a BPP, 0* is uniform random variable in the interval [0, 27r) and Ri has pdf 

fR,ir) = ^ r-in < r < rout- (39) 

Next, note that for |0j — 0o| < 2 ^ Ri ^ Rb, ^i = ^iRi has eonditional pdf 


- 


2+a|_ 

2ttu “l 


cilI^I 


_fnr < in < 

1^1 


(40) 


(a) 


When |0j—0o| > Oi in ( |40| ) is replaeed with bi, while the ol in ( |40l ) is replaeed with on when 
Rb < Ri < rout- These four eases ean be eaptured by defining a funetion V {u; [ui, 022]; c; a) as 

2+q 


75 ( 0 ;; [oil, 022 ]; c; a) = 


2tiu~ 


al^l 


{u{u — oil) — u((jJ — UJ2)} , 


(41) 


m(-) being the unit step funetion, it follows therefore that the pdf of Uj has the form 
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Using the definition in @ and the expression in ( [T7] ), we ean write the spatially averaged CCDF 
of the SINK by taking an expeetation with respeet to {U*} as 




PM=En[Pc{P,n)]=e 




i\ 


A t\ 
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t J a 


2t 




St \j=l 


To evaluate [QtX^i)], we note that the integral 
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where 


Mi (q^, c) 


Mu(x;a) = 


27rmAr(mj + ^”**^**A^/“ 

r{mi)\A\{ti\)a 

2F1 (nii + tu mi + l;mi + l + l; 


(44) 


(45) 


(46) 


{rui + f) 

rrii is as given by ( [38] ) and 2 A (a, b; c; z) is the Gauss hypergeometrie funetion. Using the 
formulation in ( |44l ) for every term in ( |42l ), E^^ [^*.(0*)] ean be evaluated so that a elosed-form 


expression for the spatially averaged CCDF of the SINK ean be eomputed from ( |43l ). Solving 
( [27] ) through numerieal integration, but with ( |43] ) in the integrand, we ean get the spatially 
averaged ergodie speetral effieieney. 


VI. Results for Random Geometry 

This seetion gives simulation and numerieal results for eoverage probability and speetral 
effieieney, whieh eonfirm the validity of assumptions made for spatial averaging. Results gener¬ 
ated using assumptions 1 to 3 are all done using a Monte Carlo simulation (whieh operates by 
randomly plaeing the interferers and bloekages aeeording to the spatial model, but then eomputing 
the eonditional outage using ([T7])). Results generated under assumption 4 ean be generated using 
either a simulation or the analytieal expression, and the methodology used will be elarified when 
the result is presented. The antenna parameters are assumed to be the same as that used earlier, 
as summarized in Table The network region A eonsidered here is an annulus with inner radius 
Tin = 0.3 m and outer radius rout = 2.1 m. The users are assumed to be randomly distributed in 
A aeeording to a BPR The simulation parameters used are summarized in Table These are 
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Fig. 9. Spatially averaged SINR coverage probability obtained via simulation for three different antenna configurations -4x4, 
16 X 1, 1 X 16 and assumption 1 with pt = 0.5. Larger Nt is advantageous and the performance is not symmetric with respect 
to Nt and Nr. 


the values used if not otherwise stated. The quantities K, W and cx^ are parameters we vary for 
eomparison later on. Varying K is equivalent to ehanging the interferer density A sinee |Al| is 
assumed to be fixed. Similarly, varying W amounts to ehanging the parameters for blockages. 
Finally, increasing captures the effect of more noise in the receiver or a lower transmit power. 

To understand the significance of scaling of the size of the antenna arrays, we plot coverage 
probability against SINR for 3 cases that have the same product of x N^. This is shown in 
Fig. 1^ where we let pt = 0.5 under only Assumption 1. As observed for the fixed geometry 


case in Section IV we see that using more transmit antennas is better than having more receive 
antennas. This is also seen in Table |V] which summarizes the ergodic spectral efficiency for 
various antenna array configurations. The reason for the asymmetrical behavior with respect to 
Nt and is that while larger Nt results in reduced probability pm oc that interferers radiate 
with larger power Gt, increasing N^ results in a decreased fraction of interferers falling within 
the receiver main-lobe which is proportional to 6^^^ oc 


We validate Assumptions 2 - 4 in Fig. 10 The plots for the CCDF of spectral efficiency with 


and without the assumptions are shown for Nt = N^ = A and Nt = N^ = IQ with pt = 1. 
We observe that, though the location of the blockages and the users are dependent in reality 
(as described by the orbital model), the assumptions of treating the blockages and users as 
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TABLE V 

Spatially averaged ergodic spectral eeficiency for various antenna configurations 


N, 

Nt 

1 

4 

16 

1 

0.6465 

1.7459 

3.2844 

4 

2.0526 

3.5963 

5.3523 

16 

3.8697 

5.5886 

7.4071 



Fig. 10. CCDF of spatially averaged ergodic spectral efficiency obtained via simulation for various transmitter and receiver 
antenna configurations Nt x N, with pi = 1. 


two independent BPPs (Assumption 2), assoeiating a distanee dependent bloekage probability 
Ph{r) (Assumption 3) and defining the LOS ball (Assumption 4) are all reasonable. The plots 


of speetral effieieney for eaeh of assumptions 1-4 when pt = l are shown in Fig. 11 

The plots in Fig. show the CCDF of the SINK obtained using the analytie expressions 


derived in Seetion V-C and eompares it with the aetual simulation results under Assumptions 1 
and 4. It is seen that the analytie expressions mateh exaetly with the setting under Assumption 
4 wherein we eonsider all the interferers within the LOS ball as unbloeked and those outside as 
bloeked from the referenee reeeiver. 

Next we look at the dependenee of the system performanee on W, the diameter of the 
bloekages. We define the throughput as the produet of pt and the ergodie speetral effieieney. The 


plot of throughput versus W in Fig. 13 shows that as W is inereased, the throughput improves. 


This is beeause, with larger W and for a fixed interferer density, the interfering signals get more 
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Fig. 11. Spatially averaged ergodic spectral efficiency from simulation when pt = 1 for various receiver antenna configurations 
and Nt = 16 with and without Assumptions 2 - 4 in Section [V-A[ 



Fig. 12. Plot showing the CCDF of spatially averaged SINR obtained from simulation and analytic closed-form expressions for 
different antenna configuration with pt = 0.7. 


blocked thus improving the SINR. The plots in Fig. 13 are for Nt = = 4 and using the 


analytie expressions derived in Seetion V-C 


In Fig.[T^ the variation of SINR eoverage probability is plotted as a funetion of A, the interferer 
density. We fix Nt = 4, N^ = 16 and pt = 1 for eomparison and use the previously derived 
analytie expressions for the plots. It is seen that as A is inereased the SINR eoverage deereases 
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Fig. 13. Spatially averaged throughput vs. W for different values of random-access prohability pt using the analytic expressions. 
Larger blockage diameter results in better throughput as the interferers are effectively blocked. 



Fig. 14. Spatially averaged SINR coverage probability vs. A for different values of SINR threshold /3 using the analytic 
expressions. Here, we let Nt = 4, Nr = 16 and pt — 1. 


rapidly initially. However, with very high density, bloeking probability also inereases, henee 
showing lower rate of deerease with inereasing A in the plots later on. This also eorroborates 
our assumption that (the users wearing) the interferers are also the souree of bloekages in the 
indoor wearables environment. 

Fig. 15 shows the variation of ergodie speetral effieieney as we vary Here we let Nt = 
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Fig. 15. Spatially averaged ergodic spectral efficiency vs for two different values of pt using the analytic expressions when 

Nt = K = 4. 


A^r = 4 and use the analytie results in Section V-C It is seen that for smaller values of the 
system is indeed interference limited as changing results in significant change in the system 
performance. However as is increased, the system becomes noise limited and different values 
for random-access probabilities of the interferers result in little change in the SINK distribution 
and hence the ergodic spectral efficiency. 


VII. Conclusion 

In this paper, we analyzed the performance of a mm Wave wearable communication network 
operating in a finite region like that inside a train car. To model the sensitivity of mm Wave 
signal propagation to the presence of human bodies in the network, we incorporated different 
path-loss and small-scale fading parameters depending on whether a link is blocked or not. It 
was seen that both interference and the probability of blockage of the interference signals are 
large when the crowd density is high, so that the SINK coverage probability decreases at a 
much lower rate with higher crowd density. We considered fixed as well as random positions 
for the interfering transmitters and assessed the impact of antenna parameters such as array 
gain and beamwidth on coverage and ergodic spectral efficiency of the system. It was seen that 
antenna main-lobe directivity and array gain play a crucial role in achieving giga-bits per second 
performance for wearable networks in a crowd. We proposed several assumptions and a model 
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to analyze the system performanee when the interferers are loeated at random loeations. These 
gave elosed-form expressions for spatially averaged eoverage probability for mmWave wearable 
eommunieation network when the user is loeated at the eenter of a dense erowd and the number 
of users is finite. The validity of the elosed-form analytie results and the assumptions were 
eonfirmed against simulations. The analytie modeling presented in this paper serves as a first 
step towards eharaeterizing SINK performanee of mmWave based ad-hoe networks in a finite 
but erowded environment, and enables one to avoid simulations to prediet performanee. 

For future work, it would be interesting to eonsider further refinements to the model ineluding 
ineorporation of 3D loeations for the deviees, and explieit modeling of refleetions of the mmWave 
signals from the boundaries of the finite network region. The work in this paper ean easily be 
extended to the ease that the referenee link ean be bloeked by the user’s own body. The proeedure 
would involve finding two eonditional outage probabilities, one eonditioned on the link not being 
bloeked by the user (using the proeedure outlined in this paper) and the other eonditioned on the 
link being bloeked by the user (adapting the proeedures so that the referenee link’s path loss is 
On and Nakagami faetor is mw). The two probabilities ean then be weighted by the probability 
of self-bloekage, whieh ean be determined based on the assumed spatial models. Using a more 
refined model to eapture this self-bloekage and ineorporating it in the analysis is an interesting 
topie for future work. 
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Appendix A 
PROOF OF LEMMAS 

The bloekages are drawn from a BPP. Consider a transmitter Xj loeated at distanee |Xj| = r 
from the referenee reeeiver. Its signal will be bloeked if there is a bloekage inside a eertain 
subregion of A, whieh we will eall the blocking region of Xj (i.e., Xj is bloeked if there is an 
objeet in its bloeking region). Sinee ^ is a eireular disk with inner radius r\n and outer radius 

if ^in < r < rout 


Tout, the bloeking region looks like Fig. 16(a) 


Y and like Fig. 16(b) 


if 


_ ry* 

out 2 ^ ^ — 'out* 
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(a) Case when r-m < r < rout — ^ (b) Case when rout- ^ <'r < rout 

Fig. 16. Figure showing the blocking region for interferer Xi with \Xi\ = r for two different cases. The blocking cone of a 
blockage Bi of diameter W that lies within the blocking region of Xi is also shown for illustration 



Fig. 17. Evaluation of area of blocking region for ri„ < r < rout — ^ 


Let Pb{i,j) be the probability that blockage Bj blocks transmitter Xi with \Xi\ = r. Since 
the blockages are placed uniformly at random, the probability that the blockage is inside the 
blocking region is equal to the ratio of the area of the corresponding blocking region and the 
overall network. For Tin < r < rout — the area of the blocking region can be evaluated as 
follows. The area can be split into regions as shown in Fig. [T^ where region Ai is a sector 
of the circle with radius rin and subtended angle = 2arcsin7^. Region A 2 corresponds 


to two identical right triangles with base length kF/2 and height w r? — (^) and region A. 


is a semicircular disk of radius fF/2. Hence, the area of the shaded region in Fig. 1 16(a) 

rW+^ 


IS 


l^il -2\A2\, where \A 2 \ = ^\jr\A - and |A| = rin^ arcsin (^) - For the 
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two cases shown in Fig. 16 we would then have, 

I ^W^-l^-\/nn^-(^)‘^-nn^arcsin(^^) + 


where 


u = 
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arcsin 


' out 


1^1 


:f)^ 


if Tin < r < rout - ^ 


if ^^out <r < rout 


(47) 


rhF 


+ rout arccos 


out 


f)V 


2\ls{s - r)(s - y )(s - ^); s = 


2rrout 

^out + r + W/2 


is the area of region Bi indieated in Fig. |16(b)[ Sinee the bloekages are independent, the 
transmitter will be bloeked if there are any bloekages — or, equivalently, will not be bloeked 
only if there are no bloekages in its bloeking region. Thus, the probability that Xi loeated at 
\Xi\ = r is bloeked is 


K 


Pb{r) = 1 - JJ(1-pb(i,j)) > 

i=i 


(48) 


resulting in the form given in Lemma 
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